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摘 要 采用 间断 有 限 元 法 (Discontinuous finite element method, DFEM) 求解 一 维 散射 性 介质 内 矢量 辐射 传输 问题 . 
推导 了 散射 性 介质 内 矢量 辐射 传输 方程 的 间断 有 限 元 离散 格式 ， 空 间 离散 采用 间断 有 限 元 离散 将 求解 域 划分 为 相互 独立 
的 单元 ,角度 离散 在 传统 的 均匀 分 段 光滑 近似 (Piecewise Constant Approximation, PCA) 角度 离散 基础 上 进行 局 部 加 密 
以 得 到 关键 方向 上 的 数值 解 . 采用 两 个 稳 态 矢量 辐射 传输 算 例 对 间断 有 限 元 法 求解 矢量 辐射 传输 方程 的 正确 性 和 计算 效 
率 进行 了 验证 , 在 此 基础 上 拓展 间断 有 限 元 法 应 用 于 求解 散射 性 介质 内 瞬 态 矢量 辐射 传输 方程 , 将 随 着 时 间 推移 过 程 中 的 
瞬 态 辐射 信息 和 最 终 的 稳 态 结果 相 比 较 验 证 了 本 文 方法 的 正确 性 . 
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Abstract The discontinuous finite nt method (DFEM) is applied to solved the polarized 


radiative transfer in a one-dimensi attering medium. The derivation in a discrete form of the 


polarized radiation governing lbns is presented, in which the spatial domain is discretized into 


finite non-overlapped disco us elements, and the angular space is discretized by the Piecewise 
Constant Approximati A) approach with a local refined modification to obtain the numerical 
results for key directions: Two steady cases are tested to verify thecaccuracy and computational 


efficiency of the developed DFEM. The DFEM is then extend t solve the transient polarized 
radiative transfer in a scattering medium, the comparisons of the@iMe-resolved Stokes vector with the 
steady ones show that the developed DFEM can handle Lo polarized radiative transfer in a 
scattering medium accurately. P X 

Key words discontinuous finite element aethon | radiative transfer; transient radiative 


transfer; discretization scheme ©: 
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F 号 R 
I Stokes 矢量 , W/(m? um sr) nax 单元 K 外 法 线 单位 向 量 No [0,27] 空间 圆周 角 离 散 个 数 
2 单位 方向 矢量 W ”权重 函数 M 离散 方向 总 数 
S 源 项 W/m? Nom 离散 单元 格 数 9 ”天 项 角 
r 空间 坐标 向 量 , m No [07] 空间 天 顶 角 离散 个 数 o BAA 
B 衰减 系数 矩阵 , m-t dg ”db6=No/r, 单位 天 顶 角 大 小 u u= cosb, RIAA HAR 
Z 散射 矩阵 Nae [x 一 d9/2, r+d9/2] 天 顶 角 局 部 加 密 离散 个 数 p BRNK 
0 引 言 


矢量 辐射 信息 相对 于 标量 辐射 信息 能 够 更 多 地 ， 程 的 求解 开拓 了 辐射 传输 理论 在 新 兴 领 域 的 广泛 应 
反应 介质 参数 特性 。 半 透明 介质 内 矢量 辐射 传输 方 M, 例如 生物 光学 工程 站 ,大气 辐射 传输 Ol, E 
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感 系统 59， 光 学 成 像 门 ， 天体 物理 学 四 等 . 近年 。 定理 进行 积分 , 为 了 简洁 省 去 表示 方向 的 上 标 n, 式 


K, 研究 学 者 们 发 展 了 若干 数值 方法 求解 半 透 明 介 (1) 可 写 为 : 
质 内 矢量 辐射 传输 , 包括 蒙特 卡 罗 法 9- 韦 ,离散 坐 
标 法 2-14) Fy 方法 M, 谱 元 法 Sl 等 . 但 这 些 方 
法 往往 具有 各 自 的 局 限 性 ， 因 此 有 必要 发 展 稳定 性 
较 好 具有 高 阶 精度 和 高 计算 效率 的 数值 方法 求解 半 Rh, none 表示 单元 K 边界 的 外 法 线 单位 向 量 . 算 
透明 介质 内 矢量 辐射 传输 ， 子 的 定义 为 : 
间断 有 限 元 法 最 初 由 Reed 和 HINI 于 1973 
年 提出 并 被 成 功 地 应 用 于 中 子 输 运 方程 的 求解 . 间 < Ag >k= |: fg dV, (fg)ox = S fgdA 44) 
断 有 限 元 在 空间 上 的 离散 单元 是 相互 分 割 的 ， 其 计 
算 区 域内 单元 之 间 的 信息 通过 边界 数值 通 量 进行 连 。 式 中 , 927 是 相 邻 单元 间 的 数值 通 量 , 本 文采 用 Lax- 
接 ， 具 有 高 阶 精度 和 较 好 的 数值 稳定 性 . 间断 有 限 。 Friedriches 数 信 通 其 格式 : 
元 方法 的 优势 使 其 受到 重视 并 开始 应 用 于 半 透 明 介 
质 内 辐射 传输 的 求解 . Cui 和 Lilis,19 率先 采用 间断 
有 限 元 方法 求解 了 半 透 明 介 质 内 标量 辐射 传输 ， 验 a I eee 示 相 邻 单 元 共同 边界 处 辐射 信息 的 
= 证 了 间断 有 限 元 方法 求解 辐射 传输 的 可 行 性 . 刘 林 值 : 
二 dpb 采用 间断 有 限 元 法 求解 了 梯度 折射 率 介 质 an 


—<I,2-VW >k +(QT -nax,W)axt 


= 3 
< BI,W >K=< S,W >k (3) 


AT = RĪ +|NEnax (5) 


Cr ”内 辐射 传输 及 半 透 明 介 质 内 瞬 态 辐射 传输 . 这 些 工 (+r), Fag -r) 6) 
= NT P 
Q 2 eh rR 边界 位 置 处 单元 K 内 部 以 及 外 部 的 辐射 信息 

,上 元 法 应 用 于 求解 半 透 明 介 质 内 矢量 辐射 作 狂 是 的 i 

E 文献 发 表 . 本 文 发 展 了 基于 Lax-Frie i A I* = lim I(rəx ~ 6nax), . 
N E cH I> = lim Froe + Sx) i 
D 射 传输 . 采用 2 PEER MITE TEARM 

A 维 散 射 性 介质 矢量 内 辐射 传 给 简 正 确 性 进行 了 检验 ， 把 式 (5) 带 入 式 (3) 可 以 得 到 矢量 辐射 传输 方 
S B e o e ns FE IA ERER: 

:二 辐射 传输 ， 将 不 同时 刻 的 辐射 信息 同 稳 态 结果 进行 | Oe 1h0-0W ont 

> TETE SER EIR AB | oo 

己 。 其 辐射 传输 的 有 效 性 mot ONT Wa + <BI,W >x= (8) 


三 ”1 辐射 传输 方程 的 间断 有 限 元 离散 i < S,W >x -5 (2. nar +|2)I7,W);x 


质 ~ 
人。 传统 有 限 元 法 需要 构造 全 局 节点 基本 孝 不 同 , 
间断 有 限 元 法 只 需要 在 每 一 离散 单元 上 构造 节点 基 
.VI(7,02")4BI(r,02")=5S(r,82") (1) 函数 . 选取 高 斯 - 切 比 雪夫 积分 节点 并 在 单元 K 上 
ccc A A 
式 中 ,了 = (I, Q, UVY 代表 Stokes 矢量 , 2 为 离散 采用 具有 5 RŽI EPR BUT, HE Stokes 矢量 


= 近似 为 : 
方向 , 6 为 衰减 系数 矩阵 ,5 为 源 项 : oe : 
S(r, 2") = BIy(r)+ T=I"(r) ~ Irol) 
M (2) = Nae = 
Z(r,Q” > T(r, R yw” BI = B(r)I”(r)~ > B IP oi(r) (9) 
n=1 ea 
对 于 间断 有 限 元 ， 由 于 其 空间 离散 节点 相互 独 oS ery DP Oe) 
WHE, HSU MER BAIL ELK = 


数 都 定义 在 每 个 基本 单元 上 . 以 任 一 单元 KK 作为 研 AF, eilr) 为 定义 在 单元 K ba eee, I? 
究 对 象 ， 选取 权 函 数 为 天 ， 在 单元 上 使 用 高 斯 散 度 。 为 第 i 个 节点 在 第 n 个 离散 方向 上 的 矢量 辐射 信 


oh Aay A iT] 
ChinaX ivq (HHF 


FH 
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E, Nex 为 单元 K 上 的 节点 数 . 将 式 (9) RAR 
(8) 中 ， 可 得 单元 K 上 的 矢量 辐射 传输 方程 的 间断 
有 限 元 离散 格式 : 


KI =H (10) 
式 中 
网 
K= [Ki]j=1, Nor;i=1, Nex) (11) 
H = [H klj=1,Nsk;k=1,4 


MEERE MAHR 豆 ”表示 为 : 


式 (10) 右边 源 项 中 的 散射 项 包含 其 他 方向 的 辐射 信 
息 ， 因 此 需要 类 似 于 离散 坐标 法 的 全 局 达 代 进行 收 
AAW. 


2 结果 及 分 析 


本 文 只 针对 一 维 问题 进行 验证 ， 以 上 公式 中 的 
方向 角 R = pi, 每 个 单元 的 节点 数 Nox = 2. 本 文 所 
有 程序 均 在 处 理 器 为 Intel(R) Core(TM) i7-4720 HQ 
CPU @ 2.6 GHz, 内 存 为 8GB 的 个 人 计算 机 上 运行 . 
2.1 散射 性 平板 介质 中 矢量 辐射 传输 

考虑 平行 光 入 射 一 维 散射 性 介质 ， 此 问题 被 称 
HAHH “L=13” H, 由 Garcia 和 Siewert!15] 提 
出 .模型 上 表面 为 透明 界面 , 下 表面 为 反射 率 p =0.1 
的 漫 反射 界面 , 界面 温度 为 OK 并 忽略 介质 发 射 . F 


K; =- < bi, R” -Voj >K + 
行 光 入 外 为 cosb. = 0.2, be = T/2， 入 射 功 率 为 
r. FRR SIS} RAK Mie 散射 性 质 的 介质 , 介 

厚度 7 = pH = 1.0, 散射 反照 率 w = 0.99, 


(2 .nor + 102°|)9i, Ojon + (12a) 


= < Bindi, Qj >K 


CD SNIDER A = 0.951， 其 散射 相 矩 阵 由 Evens 和 
O Hi =< 5,4; >K a Stephen 在 文献 [22] FAH. 

= l who) 采用 DFEM 法 对 此 问题 进行 求解 ,空间 离散 为 
LO mg (00 nox IDTE t) 20 个 均匀 单元 . 为 了 将 结果 和 和 文献 [15] 结果 进行 
©O 对 比 , 输出 量 k, Q, Us, Ve 分 别 表 示 Stokes 矢量 的 


Be w a gs SENEE 。 漫 射 量 。 首先 进行 网 格 无 关 性 验证 , 位 置 2/H = 0.5 
O oia | Ww 处 的 Stokes 矢量 L 和 V, 在 方向 余弦 p = [-1,1] 
CN —n jy aE N ; y y 


(13) 性 判断 参数 .、 图 中 给 出 了 5 种 不 同 角度 离散 条 件 下 


1A on bort < Bibi. dy> x 2/H = 04 SRD T. FI Va 的 分 布 昌 线 , 不同 角度 离散 


© sh Re we mina. 

Eo 算 子 的 定义 如 下 ， OH 1 所 示 结果 可 知 ， 当 角度 离散 个 数 达到 
二 Na \ Wc, = 40x60 BF, AAMAS BERS SCS 
ji < fis >= Df < binds > ©" 果 的 影响 不 再 显著 , 到达 网 格 无 关 性 要 求 ， 


"Ne ww 综合 考虑 计算 效率 和 精度 本 文选 取 角 度 离散 为 
Fd) = J felis $5) No x No = 40x60 计算 得 到 了 三 个 不 同 高 度 位 置 处 

a 的 Stokes 矢量 角度 分 布 ， 将 计算 结果 与 Garcia 和 
矩阵 A 可 以 改写 为 : 


=f Nsk =n 
H; = 2 S; < big1 >K- 
i=1 


Nsk 
1 a 
52" SIP (nak di, bi)art (15) 
i=l Re mi —a— Nyx N,=20x20 
1 nny WoO —e— Nyx N,=20%40 
7 |Q | ` 5 I, (Qi, bj )aK Bee ere 
i=1 —— Nyx N,=40x60 
ry ere Shee —e— Nox N,=60*100 
至 此 ， 矢 量 辐射 传输 方程 的 间断 有 限 元 离散 已 
经 完成 ， 只 需要 计算 得 到 离散 单元 的 矩阵 并 在 角度 ool 


方向 上 循环 迭代 即 可 得 到 各 节点 的 辐射 信息 。 由 于 u 
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—— N,xN, = 20x20 
—e— N,xN, = 20x40 
—a— NyxN, = 40x40 
—v— NXN, = 40x60 


È —e— NxN, = 60x100 


—0.5 0.0 0.5 1.0 
u 
图 1 网 格 无 关 性 验证 
Fig. 1 Grid independence validation 


表 1 不 同 角度 离散 条 件 下 计算 时 间 
Table 1 Computation time for different angular 
discretization schemes 


角度 离散 Ne x No 计算 时 间 Time/s 


No x Np=20x20 15.0784 

No x No=20x40 32.5988 S\ 
No x Ny=40x40 92.0667 gy 
No x Ns=40x60 174.1708 He 

No x Ns=60x100 1032.656 


Siewert 采用 Fy mares 如 图 2 
Z 


Bras. 由 图 2 可 知 本 文 计算 结 果 对 比 吻 
合 很 好 ， 验 证 了 本 文 发 展 的 疾 断 有 限 元 法 求解 散射 


性 介质 内 矢量 辐射 传输 问 
文献 [15] PRET RMA u= +0 Al u = —0 
的 辐射 信号 值 , 然而 由 图 2 可 以 看 出 采用 传统 PCA 
角度 离散 格式 所 得 结果 曲线 并 不 包含 人 = +0 和 人 
= -0 的 辐射 信息 值 . 为 了 解决 这 一 问题 , 本 文 作者 
在 传统 PCA 角度 均 勺 离散 格式 的 基础 上 对 角度 划 


REWER [x /2-d0/2, 7/2+ 40/2] ARES SIRAS V 
tt, 这 相当 于 角度 划分 总 数 为 No + Nao, 这 种 处 理 
方法 在 不 需要 大 幅 提高 计算 时 间 的 情况 下 得 到 1=0 


0.24 


0.18 DFEM This paper 
® Fy method Ref. [15] 


~ 0.12 


0.06 


0.00 
—1.0 


0.14 
ane DFEM This paper 
0.10 e Fymethod Ref. [15] 
0.08 


0.12 


DFEM This paper 
e Fy method Ref. [15] 


DFEM This paper 
® Fy method Ref. [15] 


z/H=1.0 


S 


„Č 


—0.5 0.0 0.5 1.0 


分 进行 局 部 加 密 处 理 ， 将 包含 u= 0(0 2 m/2) on 图 2 间断 有 限 元 所 得 Stokes 矢量 和 Fw 方法 所 得 结果 对 比 


Fig. 2 Comparisons of Stokes vector components by DFEM 
with those obtained by the Fy method 


附近 的 Stokes 矢量 . FX Nao = 10 对 上 述 算 例 进 行 计 
算 , 本 算 例 所 需 时 间 为 282.3374 s, 由 图 3 可 知 采用 
角度 局 部 加 密 方 法 对 角度 进行 离散 可 以 得 到 关键 方 
E y= 0 处 (辐射 信号 有 阶 跃 变化 ) 的 辐射 信息 , 这 
对 问题 求解 的 精度 具有 很 大 的 提升 ， 这 种 局 部 加 密 
角度 离散 方法 对 求解 在 某 些 特 殊 方 向 精度 要 求 很 高 
的 问题 具有 很 好 的 效果 。 
2.2 两 层 介 质 中 矢量 辐射 传输 

考虑 Evans 等 在 文献 [22] 中 所 研究 的 两 层 介质 
模型 如 图 4 所 示 ， 该 模型 上 层 是 厚度 为 4 km 的 冰 
Z, 其 衰减 系数 bi = 0.13536 km- 1， 散射 反照 率 wi 
= 0.9819; 下 层 是 厚度 为 4 km 的 雨 层 ， 其 衰减 系数 
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Br = 0.15224 km-1, 散射 反照 率 w = 0.38175, 其 他 
介质 参数 和 几何 参数 在 图 4 中 给 出 ， 


0.09 


—#— PCA uniform directions 
Nex N,=40x60 
—e— Local refine directions g 
No=40, Nag=10, N,=60 
© F method Ref. [15] f 


—s— PCA uniform directions 
Nx N,=40x60 
—e— Local refine directions 
No=40, Nao= 10, N,=60 
e Fymethod Ref. [15] 


人 < 从小 


00[ 一 -PCA uniform di oy 
NoxN, =40x60 


—0.2 -—e— Local refine directions 


Toa No=40, No=10, N,=60 
= © F method Ref. [15] 
~ 一 0.6 
一 0.8 
—1.0 
a) ‘ 
a |xQ 
0.0 0.2 0.4 0.6 0.8 N 
u 
0.2 
b : 
< —s— PCA uniform directions Ses 


NoxNo=40x60 
一 上 一 Local refine directions 
No=40, Nip=10, N,=60 
® Fymethod Ref. [15] 


0.2 0.4 0.6 0.8 


1.0 
u 


图 3 局 部 加 密 角 度 条 件 下 z/H=1 处 Stokes 矢量 分 布 


Fig. 3 Stokes vector component distributions at z/H=0 with 
local refine angular discretization 


Tyg=2.7 K 


T=245 K 
B=0.13536 km~! 

@;=0.98190 

T=273 K 


B,=0.15224 km“! 
@,=0.38175 


T=300 K 


图 4 双 层 介质 物理 模型 


Fig. 4 Physical model of a two-layer medium system 


本 文采 用 间断 有 限 元 方法 将 空间 离散 为 20 个 单 
元 , 采用 上 述 局 部 加 密 方法 进行 角度 离散 , 取 Ne = 
Nao = 10, 所 需 计 算 时 间 为 535.0904 s。 
kes 矢量 输出 为 频率 为 85.5 GHz 的 辐射 亮 
果 如 图 5 Bras. 可 以 看 出 本 文采 用 间断 有 限 元 
所 得 计算 结果 和 文献 [22] 采用 doubling-adding 算 
法 所 得 结果 吻合 得 很 好 ， 最 大 相对 误差 不 超过 2%, 
此 外 ,由 图 示 可 以 看 出 在 /= 0 处 的 Stokes RER 
有 突变 或 者 显著 的 变化 率 , H z/ 五 = 1 处 的 Q@ 在 4 
=0 会 有 负 值 出 现 . 本 文采 用 的 角度 离散 方法 可 以 很 
好 的 捕捉 到 这 些 细微 变化 趋势 和 详细 数值 信息 ， 而 
文献 [22] 中 的 数据 并 不 能 反映 出 这 些 特点 , 体现 出 
本 文采 用 计算 方法 的 先进 性 . 
2.3 瞬 态 矢 景 辐 时 传输 
考虑 辐 英 传输 方程 中 辐射 强度 随时 间 的 变化 ， 
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图 5 两 层 介质 内 矢量 辐射 传输 


Fig. 5 Vector radiative transfer in a two-layer medium system 


类 似 于 稳 态 辐射 传输 方程 ， 瞬 态 辐射 传输 方程 
在 时 刻 t=mAt 的 辐射 传输 方程 可 以 整理 为 矩阵 形 
式 : 


=n, m=n,m 


K I 


=n.m 


=H 


K, =— < Qi, R” -j >K + 
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-1) ~ apg? 
AH, Lip ”代表 离散 方向 ，m 代表 时 刻 t = 
mAt*, t* = ct/L 为 无 量 纲 时 间 , c 为 真空 中 光速 , t 
为 时 间 , L 为 模型 厚度 ,，f = 1/2 为 Crank-Nicolson 
时 间 离 散 格式 的 调节 因子 . 

采用 类 似 于 间断 有 限 元 法 求解 半 透 明 介 质 内 稳 
态 和 拓 量 辐射 传输 方程 的 计算 流程 ， 在 迭代 过 程 中 考 
虚 时 间 项 的 迭代 即 可 构造 瞬 态 矢量 辐射 传输 方程 的 
计算 流程 . 以 本 文 3.2 节 双 层 模 型 算 例 进行 验证 , 介 
质 参数 不 变 的 情况 下 界面 温度 在 t = 0 时 刻 突变 为 
图 4 中 所 示 的 界面 温度 , 介质 温度 为 0K. 采用 间断 
有 限 元 方法 进行 求解 ,空间 离散 单元 个 数 为 20， 角 


度 离散 个 数 为 Ne x No = 40x40 个 , 无 基 纲 时 间 步 
长 取 为 At = 0.1,， 最 大 无 基 纲 时 间 为 40. 

由 于 模型 为 两 层 介质 ， 分 界面 处 两 侧 的 介质 性 
质 具 有 显著 差异 , 分 析 其 投射 辐射 空间 分 布 随 时 间 
变化 , 选取 不 同时 刻 时 投射 辐射 进行 对 比如 图 6 所 
示 。 可 以 看 出 , 开始 的 一 段 时 间 内 , 介质 散射 性 质 的 
差异 使 投射 辐射 G 在 两 种 介质 交界 面 (z=4 km) 两 
侧 分 布 明 显 不 同 ， 随 着 时 间 的 推移 ， 能 量 分 布 逐渐 


趋 于 均匀 , 投射 辐射 G 的 曲线 逐渐 变 成 光滑 曲线 并 
趋 于 稳 态 结果 ,这 也 进一步 证 明了 采用 间断 有 限 元 
计 果 的 正确 性 。 


方法 求解 瞬 态 矢量 辐射 传输 结 


Incident 


图 6 投身 辐射 随时 间 变 化 的 曲线 图 


Fig. 6 “x aniation distributions at different time 


o z/H=0 Ñ z/H=1 Xb Stokes 矢量 随 
化 云图 . 由 图 示 可 知 , z/H=0 4b, Stokes 分 
7 的 值 随 着 时 间 的 推移 一 直 在 逐步 增 大 ,在 无 量 
HARFE t* = 20 左右 到 达 稳 定 ; Stokes 分 量 Q 的 值 在 
p = [0,0.5] 范围 内 出 现 了 先 增 大 后 减 小 的 情况 , Q 也 
TE t* = 20 左 右 逐 渐 趋 于 稳定 状态 . 相应 的 在 z/H=1 
Kb, Stokes 分 量 Q Æ u = [0,0.5] 范围 内 先 增 大 后 减 
小 , 但 是 幅度 没有 z/H=0 的 强烈 。 
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图 7 两 层 介 质 中 Stokes 矢量 的 时 域 分 布 


Fig. 7 Time-resolved Stokes vector components in the 


two-layer medium 


3 结 论 «~ 

将 间断 有 限 元 方法 (Discontinuous finite element 
method, DFEM) 应 用 于 一 维 散射 性 介质 内 矢量 辐 
射 传输 方程 的 求解 。 空间 离散 采用 间断 有 限 元 离 
散 ， 在 每 个 独立 单元 上 对 辐射 信息 进行 近似 ， 单 元 
之 间 通 过 包含 辐射 信息 平均 值 和 间断 值 的 数值 通 基 
连接 ,角度 离散 在 传统 的 均匀 分 段 光 滑 近 似 (Piece- 
wise Constant Approximation, PCA) 角度 离散 基础 
上 进行 局 部 加 密 以 得 到 关键 方向 上 的 精确 数值 解 ， 
给 出 了 散射 性 介质 内 矢量 辐射 传输 方程 的 间断 有 限 
元 离散 格式 . 采用 两 个 稳 态 算 例 进行 对 比 验证 ， 结 
果 表 明 本 文 发 展 的 求解 散射 性 介质 内 矢量 辐射 传输 
得 DFEM 有 很 好 的 精度 和 计算 效率 。 相 比 于 传统 
均匀 PCA 角度 离散 方法 ， 本 文采 用 的 局 部 角度 


gl 


\ 
g^ 


加 蜜 离散 方法 在 不 需要 大 幅 增加 计算 时 间 的 条 件 下 
可 以 得 到 关键 方向 的 加 密 结 果 ， 这 种 局 部 加 密 角 
度 离散 方法 对 求解 在 某 些 特殊 方向 精度 要 求 很 高 的 
问题 具有 很 好 的 适用 性 和 计算 效率 。 在 求解 稳 态 
矢量 辐射 传输 的 基础 上 本 文采 用 DFEM 求解 了 一 
维 散射 性 介质 内 瞬 态 矢量 辐射 传输 问题 ， 将 不 同 
时 刻 的 辐射 信息 和 稳 态 时 结果 对 比 验证 了 DFEM 
求解 散射 性 介质 内 瞬 态 矢量 辐射 传输 问题 的 正 
HE. 
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